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We consider a fluid queue fed by multiple On-Ofl^ flows with 
^\ • heavy-tailed (regularly varying) On periods. Under fairly mild as- 

sumptions, we prove that the workload distribution is asymptotically 
equivalent to that in a reduced system. The reduced system consists 
of a "dominant" subset of the flows, with the original service rate 
subtracted by the mean rate of the other flows. We describe how a 
(~| . dominant set may be determined from a simple knapsack formulation. 

The dominant set consists of a "minimally critical" set of On- 

Off flows with regularly varying On periods. In case the dominant 

set contains just a single On-Off flow, the exact asymptotics for the 

reduced system follow from known results. For the case of several 

On-Off flows, we exploit a powerful intuitive argument to obtain the 

►^ ' exact asymptotics. Combined with the reduced-load equivalence, the 

QQ ' results for the reduced system provide a characterization of the tail 

r~^ , of the workload distribution for a wide range of traffic scenarios. 

(<— ^ ' 1. Introduction. Over the past few decades, fluid models have gained 

■^ . strong ground as a versatile approach for analyzing burst-scale traffic be- 

^P I havior. The basic model is that of several On-Off sources, each alternating 

between activity phases (commonly referred to as bursts) and silence peri- 

j^ I ods. When active, each source generates traffic at some constant rate. 

Classical papers of Anick, Mitra and Sondhi [2] and Kosten [24] con- 
sidered a queue fed by the superposition of several homogeneous On-Off 
sources with exponentially distributed activity and silence periods. Subse- 
^ ■ quent work extended the model in various directions, such as heterogeneous 

H I source characteristics, several source states to account for various activity 

levels, or activity periods with a general Markovian structure; see, for in- 
stance, [25, 38]. Under traditional statistical assumptions, it turns out that 
the tail of the backlog distribution typically exhibits exponential decay. 
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2 B. ZWART, S. BORST AND M. MANDJES 

In recent years, empirical findings have triggered a strong interest in fluid 
models with non-Markovian activity periods. Extensive measurements in- 
dicate that bursty traffic behavior may extend over a wide range of time 
scales, manifesting itself in long-range dependence and self-similarity; see 
[26, 33]. The occurrence of these phenomena is commonly attributed to ex- 
treme variability and heavy-tailed characteristics in the underlying activity 
patterns (connection times, file sizes, scene lengths); see [5, 13, 39]. Fluid 
queues with heavy-tailed activity periods provide a natural paradigm for 
capturing these characteristics. We refer to [10] for a survey. 

Although the presence of heavy-tailed traffic characteristics is widely ac- 
knowledged, the practical implications for network performance and traffic 
engineering remain to be fully resolved. Analytical studies show potentially 
dramatic performance repercussions for infinite buffers. For moderate buffer 
sizes, though, the impact of heavy-tailed traffic characteristics is not as pro- 
nounced; see [18, 20, 30, 37]. For larger buffer sizes, flow control mechanisms 
play a critical role in preventing badly behaved traffic from overwhelming 
the buffer content; see [3]. However, the amount of backlogged traffic at the 
user, and thus the end-to-end quality-of-service, may still be significantly 
affected by heavy-tailed activity patterns. 

The effect of heavy-tailed traffic characteristics on buffer behavior also 
crucially depends on the relative amount of heavy-tailed traffic, in particu- 
lar whether or not activity of heavy-tailed flows alone can cause the buffer 
to fill. Asymptotic bounds in [15] indeed show a sharp dichotomy in the 
qualitative behavior of the workload, depending on whether the mean rate 
of the light-tailed flows plus the peak rate of the heavy-tailed flows exceeds 
the link rate or not. In case the link rate is larger, the workload distribution 
has light-tailed characteristics, whereas the link rate being smaller results 
in heavy-tailed characteristics. The exact asymptotics for the former case 
were recently obtained in [7]. For the latter case, the bounds of [15] indicate 
that one can usually identify a "dominant" set, which is a minimal set of 
flows that can cause a positive drift in the buffer. As far as bounds are con- 
cerned, all other flows can essentially be accounted for by subtracting their 
aggregate mean rate from the link rate. Somewhat related notions are de- 
veloped in [27] in the setting of M/G/co input with heterogeneous sessions. 
Exact results, however, have remained elusive for all but a few special cases. 
Results of Agrawal, Makowski and Nain [1] show that the dominance prin- 
ciple described above in fact extends to the exact asymptotics in the case 
of a single dominant flow. This may be expressed in terms of a "reduced- 
load equivalence," implying that the workload is asymptotically equivalent 
to that in a reduced system. The reduced system consists only of the dom- 
inant flow, with the link rate subtracted by the aggregate mean rate of all 
other flows. This extends results of Boxma [9] , Jelenkovic and Lazar [22] and 
Rolski, Schlegel and Schmidt [36] for multiplexing a single (intermediately) 



FLUID QUEUES WITH HEAVY TAILS 3 

regularly varying flow with several exponential flows. Related results are 
derived in [22, 35] in the context of M/G/oo input. Like the reduced-load 
equivalence, however, all these results rely on the assumption that a single 
active flow is sufficient for a positive drift in the buffer. 

In the present paper we determine the exact asymptotics for the case 
where several On-Off flows must be active for the buffer to fill (under the 
assumption that the distribution of the On periods is regularly varying [6]). 
From a practical perspective, this case appears particularly relevant, as the 
peak rate of a single flow is usually substantially smaller than the link rate. 
However, the rather subtle interaction of several flows that is involved in 
filling the buffer drastically complicates the analysis, reflecting the sharp 
demarcation in known results described above. We start with extending the 
reduced-load equivalence to the case of a reduced system consisting of several 
flows, using sample-path arguments. We then build on a qualitative under- 
standing of the large-deviations behavior to obtain the exact asymptotics 
for the reduced system. This part of the analysis is related to recent work 
of Resnick and Samorodnitsky [35] on fluid queues with M/G/oo input. 

The remainder of the paper is organized as follows. In Section 2, we 
present a detailed model description. In Section 3, we give a broad overview 
of the main results of the paper, and describe how the dominant set may 
be determined from a simple knapsack formulation. Section 4 gives some 
preliminary results. The reduced-load equivalence result is established in 
Section 5. Section 6 develops the detailed probabilistic arguments involved 
in deriving the tail asymptotics for the reduced system. In Section 7, we dis- 
cuss the relationship between the asymptotic regime considered here ( "large 
buffers") and a "many-sources" regime. 

2. Model description. We first present a detailed model description. 
We consider a queue with unit capacity (i.e., working at unit speed) fed 
by several flows indexed by the set X. For any subset E <^2, denote by 
AE{s,t) ■=J2i£E^ii^y't) the aggregate amount of traffic generated by the 
flows i £ E during the time interval (s,t]. Denote by pE := J2iGE Pi the ag- 
gregate traffic intensity of the flows i £ E (as will be specified in detail 
below). We assume p:= px <1 for stability. 

For any c > 0, £' CX, define V^(t) := supQ<j,<j{j4£;(s,t) — c{t — s)} as the 
workload at time t in a queue of capacity c fed by the flows i £ E [assuming 
V^(0) = 0]. For c> pe, let V^ be a random variable with the limiting distri- 
bution of V^(t) for t —f oo. In particular, V{t) := Vj{t) is the total workload, 
and V := Vj is a random variable with the limiting distribution of V{t) for 
t ^ oo. 

We assume the flows may be partitioned into two sets: Ti is the set of 
"light-tailed" flows; I2 is the set of "heavy-tailed" flows. For the flows i Gli 
we make the following assumption. 
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Assumption 2.1. For any c> px^, ^ > 0, 

lim x'^jV^, > x| = 0. 

The above assumption is quite weak; see, for instance, [17] for a very general 
class of arrival processes satisfying a large-deviations principle (with linear 
scaling function). However, (superpositions of) On-Off flows of which the 
activity period has a Weibull distribution satisfy Assumption 2.1 too, as can 
easily be shown using the bounds in [15] or Section 4.1 of the present paper. 
Instantaneous renewal input of which the tail of the jump sizes (bursts) is 
lighter than any power tail is covered by Assumption 2.1 as well. 

We assume the flows in I2 generate traffic according to independent On- 
Off processes, each alternating between On and Off periods. The Off periods 
of flow i are generally distributed with mean 1/Aj. The On periods Aj have a 
heavy-tailed distribution Ai{-) with mean a^ < 00. While On, flow i produces 
traffic at constant rate r^, so the mean burst size is Qirj. The fraction of time 
that flow i is On is 

_ a^ _ ^jO-i 

1/Aj + Oj 1 + AjOj 
Thus the traffic intensity of flow i is 



XiaiTi 



Pi ■=Piri 



1 + AjOj 

Before stating an important preliminary result, we first introduce some use- 
ful notation. For any two real functions /(•) and g{-), we use the nota- 
tional convention f{x) ~ g{x) to denote lim2:-»oo f{x)/g{x) = 1. Also, we use 
f{x)<g{x) to denote limsup^^^ f{x)/g{x) < 1. Similarly, f{x) >g{x) de- 
notes liminf2,.^oo f{x)/g{x) > 1. With X = Y we denote that X and Y have 
the same distribution. 

For any positive stochastic variable X with distribution function F{-), 
E{X} < 00, denote by F^'{-) the distribution function of the residual lifetime 
of X, that is, 

and denote by X^ a stochastic variable with that distribution. 

The classes of long-tailed, subexponential, regularly varying and interme- 
diately regularly varying distributions are denoted with the symbols C, S,TZ 
and ITZ, respectively (note that TZ C ITZ C S C C). Background on heavy- 
tailed distributions may be found in [16]. 

For each flow ^6X2, we assume that the On period distribution is regu- 
larly varying of index —Vi, that is, Ai{-) G TZ-Ui for some fj > 1. The next 
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result, which is due to [22], then yields the tail behavior of the workload 
distribution. 

Theorem 2.1. If A'^{-) eS, pi<c< u, then 

p{v,^ > x} ~ (1 - p,)^^f\ai > ^ j. 

c- Pi I U-c) 

3. Overview of the results. We now give a broad overview of the main 
results of the paper. As mentioned in the Introduction, asymptotic bounds 
in [15] show a sharp dichotomy in the qualitative behavior of P{V > x}, 
depending on the value of pj-^ + rjg (i.e., the mean rate of the light-tailed 
flows plus the peak rate of the heavy-tailed flows) relative to the service rate. 
In case px^ + J^Xa < 1' ^^^^ workload has light-tailed characteristics, whereas 
px-^ + rjj > 1 implies heavy-tailed characteristics. In the present paper we 
determine the exact asymptotics of P{V > x} in the latter case. For the case 
Pii +fX2<^ (where both the light-tailed and heavy-tailed input determine 
the workload asymptotics) we refer to [7]. 

3.1. Intuitive arguments. Before formulating our main theorems, we first 
provide a heuristic derivation of the tail behavior of P{V > a;}. 

Large-deviations theory suggests that, given that a "rare event" occurs, 
with overwhelming probability "it happens in the most likely way." In the 
asymptotic regime considered here ("large buffers"), the most likely way 
usually consists of a linear build-up of the workload, due to temporary in- 
stability of the system. In case of heavy-tailed distributions, the temporary 
instability typically arises from a "minimal set" of potential causes. The 
minimal set corresponds to the minimal number of causes when these are 
homogeneous in nature. In general, however, when the potential causes have 
heterogeneous characteristics, not only the number of them matters, but also 
their relative likelihood, and their relative contribution to the occurrence of 
the rare event under consideration. 

Translated to our situation, temporary instability is most likely caused by 
a "minimal set" of flows generating an extreme amount of traffic, while all 
other flows show roughly average behavior. These considerations give rise to 
the following characterization of the tail behavior of P{V > x}: 

p{v>x}~p{v^r >x}, 

with S* representing the "minimal set," and cs* := 1 — Px\s* representing 
the service rate subtracted by the aggregate traffic intensity of all other 
flows. 

We now introduce some helpful notions in order to formalize the above 
intuitive arguments. For any subset S C 2'2, define C5 := 1 — px\s as the 
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service rate subtracted by the aggregate traffic intensity of all other flows 
j £l\S. Observe that the stability condition implies ps < cs for any S CI2. 

For any subset S C T2, denote by rs := J2jes'''j ^^^ aggregate peak rate 
of the flows j G S. Define ds := rs — cs = rs + px\s — 1 as the net input rate 
(i.e., the drift) when all flows in S are On and all other flows show average 
behavior. 

A set S QI2 is called (strictly) critical if ds > (>) 0, that is, if 

rs + Pi\s> (>) 1- 

Thus, when all flows in a (strictly) critical set are On, the workload has a 
(strictly) positive drift. A critical set 5* is termed minimally critical if no 
proper subset of S is critical, that is, ds < mmj^si^j ~ Pj}- 

For any subset S C I2, denote ns '■= J2jes('^j ~ ^)- ^ strictly critical 
set S CT2 is said to be (weakly) dom,inant if /^s < (<) pu for any other 
critical set U Ql2- Observe that for a set S QI2 to be dominant, it must 
be minimally critical (because otherwise the defining property would be 
violated for any critical subset U C S). 

The quantity fis may be interpreted as a measure for the "cost" associated 
with a temporary drift ds- the probability of all flows in S being On for a 
time of the order x in steady state is roughly equal to x~^^ . Thus, a set S is 
(weakly) dominant if the flows in S being On causes the drift to be positive 
in the cheapest possible way. 

In case of light-tailed distributions, the cost minimization is usually not so 
simple; one then also needs to consider how long a certain positive drift must 
be maintained in order for a given workload level x to be reached. This issue 
does not arise in case of regularly varying On periods, since P{A[ > ax} is 
of the same order of magnitude (up to a constant) as IP{A[ > x} for any 
constant a > 1. This implies that the value of the temporary drift is not 
relevant as long as it is positive. 

Note that these heuristic arguments clearly do not hold for other subex- 
ponential distributions, such as the lognormal and Weibull distribution. In 
this case, one has F{A[ > ax} = o(F{A[ > x}), if a > 1. 

3.2. Tail behavior of the workload distribution. We now state our main 
theorem. 

Theorem 3.1 (Reduced-load equivalence). Suppose the set of flows S* C 
I2 is dominant. If Aj{-) G TZ for all j E Z2, then 

(3.1) P{V > x} ~ F{Vf: > x}, 
with 

(3.2) p{v^r >x}^(ii p,) Y. Pjoi^)^ 
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where Pjq{x) is given by (with J'i = S* \ Jq, and ds* = rs* — cs* as defined 
earlier) 



n.e^,E{A,} 



(3.3) X 



xUfI ds'Al > ^ yj{rj - Pj) + x\ JJ dyi. 
In particular, P{V > x} and Pjg{x) are regularly varying of index —fJ-s* = 

-EjGS*(2^i-l)- 

The proof of the above theorem may be found in Section 5.1 [(3.1)] and 
Section 6 [(3.2) and (3.3) and the regular variation property]. 

Note that in case the reduced system consists of just a single flow, that 
is, S* = {i*}, the tail asymptotics follow directly from Theorem 2.1. This 
is in fact the reduced-load equivalence established in [1] (under somewhat 
weaker distributional assumptions). Note that in this case the right-hand 
side of (3.2) takes the form pi* [i-*0(x) -|- Pi* {x)], with 

p,,(x)=f\ai,>^^], 
I Ti* - a* ) 

and (after a straightforward calculation) 



so that 

p,,[P^{x)+P,*{x)] = {l-p,,) P'' p(a[. > 

Ci* - Pi* I 



a* 



which is consistent with Theorem 2.1. 

In case the reduced system consists of several flows, the tail asymptotics 
cannot be obtained from known results. In fact, the analysis of the reduced 
system then poses a major challenge because of the rather subtle mechanics 
involved in reaching a large workload level. By definition, though, the re- 
duced system has the special feature that all flows must be On for the drift 
in the workload to be positive, that is, rs* — ramjQS*{fj — Pj} < cs* < rs* ■ 
In Section 6 we determine the exact asymptotics for systems satisfying this 
property, yielding the integral expression given in Theorem 3.1. 
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3.3. Knapsack formulation for determining a dominant set. We now de- 
scribe how a dominant set may be determined from a simple knapsack for- 
mulation. Recall that the On period distributions of the flows i £l2 are 
regularly varying of index — fj. 

For a strictly critical set S" C 12 to be dominant, it must necessarily solve 
the optimization problem 

sub ^rj+ J2 Pj>'^-Pii- 

j€S jei2\s 

Note that the constraint is equivalent to ds > 0. If we define 9i := ri — pi 
for all i G I21 then the above problem may be expressed in the standard 
knapsack form as 

max > [vj — 1) 

sub ^ 6j < pjj + rj2 - 1 - e, 

with U = Z2\S and e some small positive number. The above problem may 
not always have a unique solution. In case it does, the corresponding set S 
is dominant, except for the case when some set T exists which is critical but 
not strictly critical (i.e., r^ + pj\rp = 1), with p^ < P-s (see the definition 
of a dominant set). Although intriguing, this "critical case" is not further 
considered in the present paper. In this case, the temporary drift may be 
zero for a long period of time during the path to overflow. 

In case the knapsack problem has several solutions, the corresponding sets 
are weakly dominant (except for the critical case again). The next theorem 
extends the reduced-load equivalence to the case of weakly dominant sets. 

Theorem 3.2 (Generalized reduced-load equivalence; weakly dominant 
sets). Let T C 2-^2 5g ^/jg collection of all weakly dominant sets. If Aj{-) G TZ 
for all j ^ S, S" G T, then 

(3.4) P{V >x}^J2 ^{^f > ^}' 

with WiVif > x} as in (3.2), (3.3). 

3.4. Homogeneous On-Off flows. We briefly consider the case of homo- 
geneous On-Off flows as an important special case with weakly dominant 
sets. Assume that the flows i G T2 have identical characteristics. With some 
minor abuse of notation, let A[-) := j4j(-), v := fj, p := pi, r := r^, pi =p. 
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Define A^* := argminjA^ : iVr + (IJ2I - N)p > 1 - px^}. (Observe that the 
assumption px^ + rjj > 1 ensures N* < \l2\-) To exclude the critical case, 
assume that {N* — l)r + (IT2I — N* + l)p < 1 — px^ , so that the drift remains 
negative (and cannot be zero) when only A^* — 1 flows are On. 

Corollary 3.1. If A{-) g 7^, then 

P{V>x}~(^'^2j^P{V>x}, 

with 

p{v>x}~p^*5:(^*)p{i„.„.}(x), 

n=0 ^ ^ 

where -Pri,...,n}(^) ^•s given by (3.3). In particular, P{V > x} and -P{i....,n}(^) 
are regularly varying of index —N*[v — 1). 

3.5. K heterogeneous classes. We finally consider the important special 
case where each On-Off flow in X2 belongs to one of K heterogeneous classes. 
We will show how an approximate solution to the knapsack problem may 
be obtained using a simple index rule. The approximation is in fact asymp- 
totically exact in the many-sources regime. 

Specifically, consider the superposition of n On-Off flows, each belonging 
to one of K heterogeneous classes. Let a^ be the fraction of flows of class k € 
{1, . . . ,K}, with peak rate r^, mean rate pk, and an On period distribution 
which is regularly varying of index — ffc. Let the service rate be n (instead 
of 1), and let V*^") be the stationary workload. The knapsack problem then 
takes the form 



K 



mm 



'^nk{vk 



nfeG{0,...,nafe}^^_l^ 

K K 

sub ^ nfcrfc + ^ {nau - nk)pk > n. 
fc=i fc=i 

Unfortunately, the above problem cannot be easily solved due to the inte- 
grality constraints. Intuitively, however, one may expect that as n grows 
large, the integrality constraints should have a negligible effect, so that a 
continuous relaxation with Uk G [0, nak] should give a good approximate 
solution. 

This relaxation may be solved using a simple index rule. Index the K classes 
in nondecreasing order of the ratios 

Ik ■■= i^k - l)/(?'fc - Pk)- 
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For any k £ {1, . . . , K}, define cjfc := Em=\ (^mrm + Ylm=k o-mPm- Determine 
the (unique) index / such that 1 € {ai^i,a]\. Then take n^ = na^ for all 
classes k <l, n1 = for all classes k > I, and n* = n(l — ai^i)/{ri — pi). 
This yields the (crude) approximation 

(3.5) P{V(") > x} « x""^, 

with /i := X)/—! '^fc(^fc ~ 1) + (1 ~ <yi~i)li- In Section 7 we prove that the 
above approximation is logarithmically exact in the many-sources regime. 
In particular, one may show that the limits for x ^ oo and n — > oo commute 
if one considers logarithmic asymptotics. 

Theorem 3.3 (Robustness of logarithmic asymptotics). 

1 logP{V(") >nx} , , 1 logPlVW > nx] 
lim lim = lim lim . 

n^oax^co n lOgX x^oon^ooTT, logX 

The proof of the above theorem may be found in Section 7. Although log- 
arithmically exact, the approximation (3.5) may not be appropriate from a 
practical perspective. In particular, it is shown in Section 7 that an ana- 
logue of Theorem 3.3 cannot hold if one considers exact asymptotics. This 
"negative" result is reminiscent of a phenomenon occurring in heavy-traffic 
theory where two limiting regimes lead either to stable Levy motion or to 
fractional Brownian motion; see, for example, [31] and references therein. 

4. Preliminary results. In this section we collect some preliminary re- 
sults which will be used in later sections. 

4.1. Bounds. We first derive some simple bounds for the workload dis- 
tribution P{V^ > x} for subsets S" C T2. 
For any subset S CI2, c<rs, define 



„v-c ^ 



X 



The next lemma gives a lower bound for P{V|. > x} which may also be 
found in [12]. 

Lemma 4.1. Let S <^l2- For c<rs, 

F{V's>x}>P§{x). 

Proof. Consider the event that at some arbitrary time t all flows j G 
S have been On since time t ^— or longer. This event occurs with 



FLUID QUEUES WITH HEAVY TAILS 11 

probability Pg{x), and implies that the workload at time t is larger than 

rg—c rs — c 



For any subset S CI2, c<rs, define 



^s:=n 



r, 



Pj 



Tj -pj+c-rs 



The next lemma establishes an asymptotic upper bound for P{V^ > x} for 
the case where S is a minimally critical set with respect to the capacity c. 

Lemma 4.2. Let S (^12- Ifce{rs-m.mj(zs{rj-pj},rs), and A'^j{-) eS 
for all j € S, then 

F{V's>x}<K^sPsi^). 

Proof. For any i G S, denote di := c — rs + ri. Observe that di > pi 
since c> rs — {vi — pi). We apply the usual technique to obtain an upper 
bound: split the capacity. Formally, we have the sample-path upper bound 

(4.1) Vl{t) < Vf^{t) + y;(J«(t) = Vt^{t) 

for all i € 5. 

In the stationary regime, using Theorem 2.1, 

¥{V% >x}< ¥{\f > X for all j G S} 



: n P{V;^ > X} 

ie5 



.n'-«>^'H>^} 



= n pj — ^-^-^ — H A^' > ^— 

j<fs ^j- Pj+c-'^S I ■' rs-c 
= K-sPl{x). D 

Corollary 4.1. Let S (^12- If c £ {rs - minjfz s{rj - pj},rs), and 
A'^j{-)£SforalljeS, then 

Plix)<¥{V's>x}<K^sPsi^). 

Proof. The proof follows directly by combining Lemmas 4.1 and 4.2. 
D 
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Corollary 4.2. Let S CI2. If A'j{-) G ITZ for all j G 5, then, for any 
closed interval T C (^rs — Taiiij^sWj ~ Pj}jfs)^ there exist constants K^^' , 
K^"^' independent of c, such that for all cGT, 



K^Ps^x) < P{V^ >x}< K^^^Psix) 



with 



Ps{x):=Y[F{A'j>x}. 

Proof. The statement follows directly from Corollary 4.1 and the fact 
that A'^ji-) e in CS for ah j eS when observing that A'j^-) e ITZ, j e S 
implies that 

hmsup % <cx), 

x— >oo J^ (X) 

if ci,C2 G T. D 

We now derive some general bounds for the total workload distribution 
P{V > x} which will be crucial in establishing the reduced-load equivalence. 

For any c > 0, E (ZI, define Z'^{t) := supQ<j,<j{c(t — s) — A^(s,i)}. 
For c < pe, let Z^ be a random variable with the limiting distribution of 
Zj^{t) for t ^ 00. Let Q C 2^ be the collection of all minimally critical sets. 

We first present a lower bound. The idea is as follows: V^ being large 
for some minimally critical set i? G A basically implies that V must be 
large too, unless the other flows j ^ E persist in below-average behavior. 
Excluding such below-average behavior (reflected in large values of Zj> ^) 
from the event {V > x} yields the following lower bound for P{V > x}. 

Lemma 4.3. Let A C il. Then for any 6>0 and y > 0, 

P{V >x}>Y^ ¥{Y%^+' >x + y}P{Zf;^"~' < y} 



EeA 



E n ^\^r >- 

Ei,E2eA,EiytE2 J&E1UE2 



Proof. Sample-path wise, 

V{t) = sup {A{s,t) -{t- s)} 

0<s<t 



: sup {AE{s,t) + Aj\E{s,t) - {CE + 6){t - s) 
0<s<t 

-ipj\E-m-s)} 
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> sup {AE{s,t)-{cE + 6) (t-s)} 

0<s<t 

+,HM^\Eis,t) - [px\E -m-s)} 

0<s<t 

= sup {AEis,t)- {CE + S){t- s)} 

0<s<t 

- sup {{px\E - S){t - S) - Aj\E{s,t)} 
0<s<t 

for all E gA. 

In the stationary regime, for any 6 > and y > 0, using the independence 

of V'i^' and Zf^\f -', 

P{V > x} 

> P{Vf ^-^ - Z^^^^^^"^ > X for some ^ G A} 

> F{V£^^ >x + y, 7f^^E~^ < y for some E e A} 
(4.2) > P{V^^+'^ >x + y, t!'-^\e~^ < y for exactly one E e A} 

- E mT'>-+y,^?\'Er'^y^ 

>Envr'>^+2/}nz^^r'<y} 

- E mT'>^,^T'>^}. 

Ei,E2i=A,Ei¥=E2 



As in (4.1), 
(4.3) 






for all i^E. 

Note that ce — ''^E\{i} > Pi foi^ all i G £', E ^ A, since E is minimally 
critical. 
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Hence, 



y^^+'^(t) < yf +'^(t) 



for allieE, E € A. 
Thus, 



F{V''1^'-'>X,V2^-''>X} 



(4.4) 



El ^ ^1 * E2 

< PlVj^"^*^ > X for all j G ^i, VJ*^"^^ > x for all j G £^2} 



h<5 



P{V7"^ > X for all jeEiUE2} 



J 



= n ip{vf+ >x}. 

Substituting (4.4) into (4.2) completes the proof. D 

We now provide a corresponding upper bound, which is somewhat more 
involved. The idea is as follows: V being large essentially means that V^ 
must be large for some minimally critical set i? G A too, unless the other 
flows j ^ E exhibit above-average behavior. Extending the event {V > x} 
with possible above-average behavior of the flows j ^ E (manifesting itself in 

large values of V^-f^ ) leads to the following upper bound for P{V > x}. 

Lemma 4.4. Let A C il. Then for any 6,e > sufficiently small and y, 

P{V >x}<Y, IP{V^^~^ >x- y} + F{Vjl^^' > x/AT} 
EeA 

+ E n^i\E^' >y}I[ n^?^' > x/M} 

EeA jeE 

+ E nnvf^%x/AA}, 

EGfl\AjGE 

with J\f := \T\ denoting the total number of flows. 

Proof. As before, we divide the capacity to obtain the sample-path 
upper bound 



V{t)<V^^-~\t) + V^^f^\t) 



'I\E 
for all £■ G A. 



In addition, for e > sufficiently small, V{t) > x implies V^^ ^ (t) > x/J\f, 

rP. 
3 



or there exists a minimally critical set S G such that K- ^ (t) > x/J\f for 



all j G S. 
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This may be seen as follows: suppose that it were not the case, that 
is, Vj '^ (0 — x/J\f, and for every minimally critical set S" € ^2 there ex- 
ists a j (depending on S) such that V^'^ ^{t) < xjM . Then the set Jit) := 

{j £ ^2 ■ Vj ^ (t) > x/J\f} does not contain any minimally critical set, hence 
rj{t) + Pi\j{t) < 1- This means that Pi\j{t) +Me<l- rji^t) for e > suffi- 
ciently small. Thus, noting that Pi\j{t) = Px^ + Px-i\j(t)-, 






jGX2\J(i) 

<\X\J{t)\xlM 

< X, 

contradicting the initial supposition. 

In the stationary regime, for any 6,e > sufficiently small and y, using 
independence 

P{V > x} 

< P{Vf -^ + V^^^^"^"^'' > X for all E e A, 

^S'^' > x/M or V^^+^ > x/M for all j G 5 for some 5 G f]} 

< F{V^^"^ > X - y or V^^^^"'^^ > y for all ^ G A, 

Vj^'^^ > ^1^ or Vj^^^ > x/M for all j G 5 for some 5 G f^} 

< Y. P{V|^-^ >x-y}+ P{V^^i+' > x/N] 
EeA 

+ Y, P{Vj^^+^ > x/M for all j G S, ^^^e^^ > y for all EeA} 

Sen 

< E P{Vf "^ >x-y}+ P{Vj^^+^ > x/M} 

EGA 

+ E F{Vf +" > x/J\f for ah j G E, V^^^f "^^ > y} 
EeA 

+ E IP{Vj'^^>x/A/'for aU jG^} 

£60\A 



16 B. ZWART, S. BORST AND M. MANDJES 

< Y. P{V|^-^ >x-y}+ P{Vj^+' > x/A/"} 
EeA 

+ E n^T\'E^' >y}I[ ^v^?^' > ^/^} 

EeQ\AjeE □ 

4.2. Stationary workload representation. We now give a convenient rep- 
resentation for the stationary workload V|; , with -E C T2 an arbitrary set of 
heavy-tailed On-Off flows. We start from the definition V^{t) := supo<s<t{^£;(si t) ■ 
c{t — s)} [assuming V^(0) = 0]. Since the process Ae{-, •) has stationary and 
reversible increments, we have 

sup {AE{s,t) — c(t — s)} = sup {Ae{0, s) — cs}. 

0<s<t 0<s<t 

In the sequel, we simply use the latter expression as the definition of V^(i). 
Accordingly, for c> pE, the stationary workload as t ^ 00 may be repre- 
sented as 

V|:=sup{^£(0,t)-ct}. 
t>o 

Explicit constructions of Ai{0,t) (satisfying the stationarity condition) may 
be found in [15, 19]. For completeness, we review the construction in [19], 
which will be extensively used in Section 6. 

Let {Aim, m > 0} be a sequence of i.i.d. random variables representing On 
periods of flow i. Similarly, let {Uj^rn""^ ^ 1} be Off periods. Define three 
additional random variables A[q, U[q and Bj such that A[q = A[, U^q = U[ 
and 

Note that B^ = 1 corresponds to flow i being On (in stationarity). 

To obtain a stationary alternating renewal process, we define the delay 
random variable Djo by 

D,o = B,A[o + (1 - Bi){Vlo + Aio). 

Then the delayed renewal sequence 

{Zin,n> 0} = J Dio,Dio + Y. (Uim + Ai^),n > 1 i 

I m=l J 

is stationary. 
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Next, we define the process {Ji{t),t > 0} as follows. Jj(t) is the indicator 
of the event that flow i is On at time t. Formally, we have 

Ji{t) = Bi]l{t<A-j + (1 - Bi)l{U[o<i<u-o+A,o} 

oo 

+ 2^ l{Zi„+Ui,„+i<t<Zi,„+i}- 
n.=0 

The On-Off process {Ji{t),t > 0} is strictly stationary; see [19], Theorem 
2.1. The process {Ai{0,t),t > 0} is defined by 

Ai{0,t) :=ri / Ji{u)du. 
Jo 

Finally, note that the number of elapsed Off periods during [0,t] which 
started after time is given by 

(4.5) iV^(t) := max{n : Zi,„_i + Ui„ < t}. 

We conclude this section with the following useful lemma. 

Lemma 4.5. LetSC'l2- IfAj{-) G TZ for all j e S andc£ (rs-minjgsjrj 
Pj},rs), then 

,. ,. P{supt>MjAs(0,t)-(c-e)t}>x} 
hm iimsup = -— r = U, 

for any e S [0, rs — c). 

Proof. For t > Mx, write 
As{0,t)-{c-e)t 

= As{0, Mx) - (c - e)Mx + As{Mx, t)-{c- e){t - Mx), 

and observe that As{Mx,t) = ^^(O,* — Mx) since the process ^^(O,*) is 
stationary. Thus, for 6 > sufficiently small, 

p( sup {As{0, t)-{c- e)t} > x\ 

"l sup {As{0, Mx) - (c - e)Mx 

+ As{Mx, t)-{c- e){t - Mx)} > x\ 
As{0,Mx)-{c-e)Mx 



{t>Mx 
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+ sup {As{Mx,t)-{c-e){t-Mx)}>x\ 

t>Mx } 

< ¥{As{0, Mx) - (c - e)Mx > -6{c - e)Mx} 

+ p( sup {As{0, t - Mx) - (c - e)(i - Mx)} > (1 + 6{c - e)M)x\ 

= F{As{0, Mx) > (1 - (5)(c - e)Mx} 

+ p( sup {As{0, t - Mx) - (c - e)(t - Mx)} > (1 + 6{c - e)M)x\ 

< p(sup{^5(0, t) - (1 - 25){c - e)t} > 6{c - e)Mx\ 

U>o J 

+ p(sup{A5(0, t) - (c - e)t} > (1 + 5{c - e)M)x\ 
U>o J 

= p{v(|-2^)(^-=) > 6{c - e)Mx} + P{V^--' > (1 + 6{c - e)M)x}. 

Using Corollary 4.2, for 5 > sufficiently small, 

Hsupt^Mxi^sjO, t)-{c- e)t} > x} 
P{V^ > x} 

K^^^Ps{6{c - e)Mx) K^'^'^Ps{{l + 6{c - e)M)x) 



K^^)Ps{x) K(^)Ps{x) 

Now let X — > oo and then M ^ oo [use the fact that Psi-) is of regular varia- 
tion] . 

D 

5. Reduced-load equivalence. In this section we provide the proofs of the 
various reduced- load equivalence results stated in Section 3. The proofs of 
the complementing results for the reduced system are presented in Section 6. 
Section 5.1 treats the case of several weakly dominant sets, culminating in a 
proof of (3.4). See Theorem 5.1; this also gives the proof of the special case 
in Theorem 3.1. In Section 5.2 we extend the results to the case of additional 
instantaneous, heavy-tailed input. 

5.1. Proof of (3.4). Recall that T denotes the collection of all weakly 
dominant sets, and that f2 represents the collection of all minimally critical 
sets. 

Assumption 5.1. For any y and 5 > 0, 

•^^ ^ ^->oo P{V^ > x} 
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is independent of y. In addition, limsio Fg{5) = 1. 
Assumption 5.2. For any y and 5 > 0, 

G-,W:=hmsup ^f^. ^ ^^ 
is independent of y. In addition, lim^^o ^'5(5) = 1. 
Assumption 5.3. For any e > 0, 

X^^^D P{V^ > X} 

Assumption 5.4. For any e > 0, 

Hsie) := hmsup ^^^^^-^^ < 00. 

Assumption 5.5. For any pair of sets S £T, E efl\T, ior any e > 0, 
lim — ^ ^ = 0. 

x^oo P{V^ > x} 

Theorem 5.1 (Generalized reduced-load equivalence; weakly dominant 
sets). Suppose the sets S gA satisfy Assumptions 5.1-5.5. Then 

P{V >x}'^J2 PV^f > ^}- 

sga 

Proof. As before, the proof consists of a lower bound and an upper 
bound which asymptotically coincide. For compactness, denote Q{x) := 

(Lower bound.) From Lemma 4.3, for any 6 > and y > 0, 

P{V > x} > ;^ P{V^^'+^ >x + y}P{Z^^\f '^ < y} 

5GA 

- E n n^'/'-'yx/Af}. 

Si,S2eA,Sij^S2jeSiUS2 

Note that if 81,82 G A, 81 7^ 5*2, then 81 U 82 cannot be a minimally critical 
set, so that ^i U 5*2 ^ A. 

Thus, using Assumptions 5.1 and 5.4, and the inequality 

— -— > mm — 

l^i Oi I hi 
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for ai,bi > 0, we obtain 

P{V > x} 



liminf- 
x^oo Q{x) 



T-CS+S 



P{V^'5+'' >x + y} 



x— >oo 

seA 



> liminf y P{Z^f.^ < y} — - 

2_] li™ sup 



n,e5.u5.]P{Vf^^>^/A/'} 



> liminf minP{Z^f,^-' < y}EiXL^^±ll 






^/...„.:„.nvr'>^+2/} 



> min P|Z!;:^j' < y} lim inf '^ '= ^, ^ 

= mmF^-(<^)P{zf^^/-'<y}. 

Letting y — > oo, then 5 J, 0, we obtain 

, PjV > xj 
liminf ^ / >1, 

which completes the proof of the lower bound. 

(Upper bound.) From Lemma 4.4, for any (5 > and y, 

P{V > x} < y P{V^«-'' >x-y} + P{V^^i^' > x/A/"} 
5eA 



+ E ip{v^^^/^' > y} n n^7^' > x/AA} 



5eA iG5 

EGii\AjeE 

Thus, using Assumptions 5.2-5.5, and the inequality 

— -— < max — 

22i Oi ■>- bi 

for ai,bi> 0, 



P{V > x} < hmsup ^ ^ ',. ^ + hmsup ^^ 

x^oo (J , h;l,''^J X— >oo V(,xj 

+ E ip{ Vj\s > y} 1™ sup -^-- 
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, ^ ,. niei?F{Vf+%x/A/-} 

< lim sup max ■ 



X- 



r^^sGA p{v|s>x} 






SeA 

-5 



^ S'-T^^fe^ + S ^s(.)P{v-r' > .} 



s ^-^J 5eA 



maxG^^(5) + Y: ^5(e)P{vJ^'/''' > y}- 



SeA 
Letting y — > 00, then (5 | 0, we obtain 

limsup — -— - — < 1, 

which completes the proof. D 

In order to complete the proof of the reduced- load equivalence result (3.1), 
it remains to be shown that a dominant set 5* C X2 with Aj{-) G TZ for 
all j € S* satisfies Assumptions 5.1-5.5. That is done in the following two 
propositions for 5 = 5*. 

Proposition 5.1. Let S CI2. If Aj{-) e TZ for all j e 5, then Assump- 
tions 5.1 and 5.2 are satisfied for any c G {rs — niinj^s{fj — Pj}i^s)- 

Proof. We first prove that Assumption 5.2 is satisfied. It follows from 
Theorem 6.4 (see also Corollary 6.1; it is important to note here that results 
from Section 6 do not rely on the results of this section) that if Aj(-) € TZ 
for all j G 5, then P{V^ > x} G TTZ. Since ZTZ C £, it suffices to prove that 
the assumption is satisfied for y = 0. 

Let £ G [0,rs - c), and let 5 G (0,e]. Then 

H^s'^ >x} = p|sup{A5(0, t)-{c- 5)t} > x] 
U>o J 

<p( sup {As{0,t)-{c-S)t}>x\ 
+ p( sup {Asi0,t)-{c-5)t}>x\ 
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<p( sup {As{0,t) - ct} > {1 - 6^/^)x\ 
+ p( sup {As{0,t)-{c-e)t}>x\. 



Thus, 

P{V".~^>a;} , ¥{V% > (1 - 6^/'^)x} 
lim sup —^-^ r^ < lim sup — ^^ — ^—-i , 

IP{suPi>,5-i/2{^s(0, t) - (c - e)t} > x} 

+ lim sup = ;;rrrT7 ^ • 

x^oo^ P{V^ > x} 

The fact that PjVI' > x} G T7^ imphes that the first term tends to 1 as 5 J. 0, 
while Lemma 4.5 (with M = S~^''^) shows that the second term then goes 
toO. 

The proof that Assumption 5.1 holds is similar, and therefore omitted. 
D 

Proposition 5.2. Let S Ql2- If Aj{-) e TZ for all j G S, then Assump- 
tions 5.3 and 5.4 are satisfied for any c> ps- If, in addition, S is a weakly 
dominant set, then Assumption 5.5 is satisfied as well. 



Proof. Using Lemma 4.1, 

p{v^>x}> n^jip|Ai> 



Assumption 5.3 then follows from combining Assumption 2.1 and the as- 
sumption that Aj{-) G TZ for all j G S. 
Theorem 2.1 gives 

P{vf +" > x/M} ~ (1 -p^)PiAa] > ^•^_ 

e [ rj pj e 

for all j G ^2. 

Assumption 5.4 then follows from the assumption that Aj{-) G TZ for all 
J G S, and so does Assumption 5.5 in case 5 is a weakly dominant set. D 

5.2. Additional instantaneous input. So far we have considered a sce- 
nario with only fluid heavy-tailed input. We now extend the reduced- load 
equivalence to the case with additional instantaneous, heavy-tailed input. 
We thus allow for an additional subset of flows I-^ (^I which generate in- 
stantaneous traffic bursts according to independent renewal processes. The 
interarrival times between bursts of flow i are generally distributed with 
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mean l/Aj. The burst sizes Bj have a heavy-tailed distribution Bi{-) with 
mean /3j < oo. Thus the traffic intensity of flow i is pi := XiPt. 

For each flow i G X3 , we assume that the burst size distribution is regularly 
varying of index — fj, that is, Bi{-) G 'R—ui for some Vi > 1. The next result 
which is due to [32] then gives the tail behavior of the workload distribution 
for a single flow i € X3 served in isolation. 

Theorem 5.2. // -B[(-) eS, pi< c, then 

(5.1) P{Vf > x} ~ ^^P{B[ > x}. 

C Pi 

In order to formulate the results, we need to extend the concept of dominance 
introduced in Section 3. A flow i GI3 is said to (weakly) dominate a flow 
j €l3 if t'i < (<) i^j. A flow i G X3 is said to (weakly) dominate a critical set 
S CI2 ii Ui — 1 < (<) J2j£si^j — !)• A critical set S C X2 is said to (weakly) 
dominate a flow i € X3 if i/j — 1 > (>) J2jesi'^j ~ !)• 

A flow i € X3 is called (weakly) dominant if it (weakly) dominates all other 
flows j S 23 as well as all critical sets /S C X2. A critical set 5^X2 is called 
(weakly) dominant if it (weakly) dominates any other critical set C/ C X2 as 
well as all flows j E X3 . 

Theorem 5.3. Let /C C X3 and T C 2^^ be the collection of all weakly 
dominant flows and all weakly dominant sets, respectively. If Bi{-) E 7^ for 
all i E /C, and Aj{-) G TZ for all j G S, S* G T, then 

(5.2) P{V > x} ~ ^ F{V^'' >x}+J2 ^{^7 > ^}' 

ie/C 5GT 

with P{V^' > x} and P{V^® > x} as in (5.1) and (3.2), (3.3), respectively. 

The proof of the above theorem is similar to that of Theorem 5.1 after a 
few modifications to Lemmas 4.3 and 4.4. 

It may be worth mentioning that Theorem 5.3 continues to hold under 
the condition Bl{-) G S for alH G /C, provided there are no weakly dominant 
sets of On-Off flows (the concept of dominance may be extended to subex- 
ponential distributions in a straightforward way). In particular, when there 
are simply no On-Off flows at all, one obtains the extension of Theorem 5.2 
to the single-server queue fed by a superposition of renewal processes (which 
is not a renewal process). This result was obtained as Theorem 4.1 in [4], 
using a different approach. 

Theorem 5.3 also provides an extension of a recent result in [11], who 
study an M/G/1 queue with two different speeds of service using an analytic 
approach. A queue with two service speeds fits into our framework by the 
observation that the varying service capacity can be regulated by an On-Off 
source. 
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6. Tail asymptotics for the reduced system. In this section we derive 
the tail asymptotics for the reduced system. In particular, we give a proof 
of (3.2) and (3.3). 

For notational convenience, let c be the capacity of the reduced system, 
let the set of flows be indexed as ^7 = {1, . . . ,N}, and denote r := rj and 
A{0,t) := Aj{0,t). By definition, the reduced system satisfies the following 
two properties: 

(i) The On period distribution of flow i is regularly varying of index — fj < — 1, 
that is, Ai{-) £ Tl-u,- 

(ii) All flows must be On for the drift of the workload process to be 
positive, that is, c£{r — minj=i^...^Ar{rj — pi}, r). 

We now state our main theorem. 

Theorem 6.1. Consider a queue of capacity c fed by N On-Off flows. 
If c £ {r — minj=i^.,.^7v{rj — Pi},r) with r = X]i=i^ii ""■^ ^j(') ^ ^ /^^^ ^^^ 
J = 1, . . . , A^, then 



V=l / JoC{l,...,Ar} 



{1,...,7V} 

where Pj^^ (x) is given by (with Ji = {1, . . . , N} \ Jq) 
(6.1) X 



/ n ^1 (^ - c)Ai > E yi(^i - Pj) - (^ - c)yi + X \ 

J2/,6(0,oo),iGJi .g^^ I jeJi J 

X n ip| (^ - c)Ar > E %(^j - pj) + ^ [ n dyi- 

An asymptotic characterization of Pjq{x) which may be useful for further 
analysis is provided in Section 6.4. This characterization also shows that 
P{ V'^ > x} and Pj^ [x) are regularly varying, and gives an expression for the 
pre- factor in the asymptotic expansion of PjV^ > x}. 

The remainder of this section is organized as follows. Detailed heuristic 
arguments are given in Section 6.1. In Section 6.2, we prove some preliminary 
results on the most probable behavior of the process {^(0, t) — ct}. The proof 
of Theorem 6.1 is then completed in Section 6.3. Section 6.4 deals with the 
asymptotic behavior of Pjq{x). 
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6.1. Heuristic arguments. The proof of Theorem 6.1 is quite lengthy. 
Nevertheless, it is based on a simple intuitive argument: the most likely way 
for V^ = supj>Q{A(0,t) — ct} to reach a large value is that all flows have 
been simultaneously On for a long time. Specifically, each flow is likely to 
contribute through exactly one "long" On period; apart from these long On 
periods, the flows show typical behavior. 

The above heuristic argument may be used for computing sup(>o{^(0, t) — 
at}. Let us say that the long On period of flow i begins at time Sj and ends 
at time Sj + 1 j . Define 

t* := min {sj + ij} 

i=l,...,N 

as the time epoch at which the first of the long On periods finishes. One 
may also interpret t* as the time epoch at which the process {A{0,t) — ct} 
reaches its largest value. Note that Aj(0, Sj) ss p^Si, Ai{si,Si + tj) = rjtj, and 
Ai{si + ti,Si + tj + t) ~ Pit, t >0. One thus obtains, using the fact that 
c€{r- mini=i^...,Ar{ri - pj, r), 

sup{A(0, t) - ct} PS A(0, t*) - ct* 

N 

(6.2) f«^[piSi + ri(t*-Si)]-ct* 

4=1 

JV 



^{Pi - ri)si + (r - c)t* 



i=l 

The problem is thus reduced to calculating 




(6.3) F{}^{pi-ri)si + {r-c) mm {si+ti} > x} . 

t=l,...,N J 

Although the proof is based on the representation V^ = sup^>g{^(0,t) — 
ct}, it is useful to keep the original workload process supo<s<t{^(s,t) — 
c{t — s)} in mind as well. Figure 1 shows a typical scenario leading to a large 
workload level (so small fluctuations are ignored) in the case of two On-Off 
flows. 

At a certain time ujq, the first long On period begins. Before that time, 
both flows show average behavior. The queue starts to build (at rate ri + 
r2 — c) at time wi when the second long On period begins, and reaches its 
largest level at time co^. Level x is crossed at time uj2- 

Between times ^3 and ll'4, the queue drains at rate c — ri — p2- flow 1 is 
still in the middle of its long On period, and flow 2 shows average behavior 
(remember small fluctuations are neglected). The process is still above level x 
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V-(f) 




Fig. 1. Typical overflow scenario for two On~Off flows. 



between times u;4 and ws. However, here both flows show average behavior 
again, causing a negative drift c — pi — p2- 

The figure illustrates why the analysis of the reduced system is still quite 
complicated: 

(i) Although the long On periods must significantly overlap, the differ- 
ence between the finishing times of these On periods can be quite large (of 
order x, hence not negligible). 

(ii) Given that the observed workload is larger than x, it is not necessarily 
the case that all flows are in the middle of their long On periods. In Figure 1, 
this is only the case in the time interval (a;2,W3). In fact, for any given flow, 
its long On period may have finished a long time ago. Consequently, there 
are 2 different possibilities (corresponding to which subsets of the flows 
are still in the middle of their long On periods). Sample-path wise, there are 
A^ + 1 different time intervals in which the workload may be larger than x 
(depending on how many of the flows are still in the middle of their long On 
periods). 

(iii) Specifically, given that the observed workload is larger than x, it may 
still have been even larger at an earlier time epoch. In Figure 1, this is the 
case in the time intervals (u;3,ti;4) and (ti;4,a;5). 

These complications do not arise if one considers a related problem, which 
concerns the overflow probability in a fluid queue with a finite buffer of 
size X. As is shown in [23], the analysis of the reduced system is then consid- 
erably simpler. It suffices to use bounds which are similar to Lemmas 4.1 and 
4.2, and to combine these with the asymptotic results for a single On-Off 
flow in [21, 40]. 
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6.2. Characterization of most probable behavior. In this section we prove 
some preliminary results characterizing the most probable behavior of the 
process {A(0,t) — ct} given that it reaches a large value. In particular, we 
formalize the following two heuristic statements, resulting in a formal version 
of (6.2). 

(i) Each flow contributes to sup(>o{yl(0,t) — ct} through exactly one 
"long" On period. 

(ii) Apart from these long On periods, the flows show typical behavior. 

An On period is referred to as "long" when larger than ex, with e some 
small but positive constant. In order to formalize the above statements, we 
need to keep track of how many such long On periods occur. 

With that in mind, we define Afi{A, B), for intervals A,B CI [0, cxd), as the 
number of On periods of flow i of which the length is contained in A and 
which overlap (in time) with B. For compactness, denote 

Af^{u,t)=Mi{{u,oo),[0,t]). 

We now proceed with a few preparatory lemmas. 

First we show how to obtain an upper bound for the workload process 
in terms of a simple random walk. As in (4.1), we have V^{t) < V^ '(i) 

for all i = 1,. . . ,N, with di := c — rjxux = c — r + ri. Recall that V- '■{t) = 
supo<s<t{^i(0,s) — dis}. Now let, for fixed i, S^ = Xji + ■ • ■ + Xj„ be a 
random walk with step sizes Xj^ = {^i — di)Aim — diVim, with Aim and 
Ujm i.i.d. random variables distributed as the On and Off periods of flow i, 
respectively. 

Since cG {r — minj=i^...^Ar{rj — pi}, r), we have pi < di for alH = 1, . . . , A^, 
so that E{Xji} < 0, that is, the random walk has negative drift. Because of 
the sawtooth nature of the process ^^(0, s) — diS, we have 

sup {Ai{0,s) - dis} < {ri - (ii)(BjA^o + (1 " 'Bi)Mo) + sup Sin, 

0<s<t n<N^{t) 

with N^{t) denoting the number of Off periods of flow i elapsed during [0, t] 
which started after time [for a formal definition see (4.5)]. 

The above observations are summarized in the following auxiliary lemma. 



Lemma 6.1. For all e> 0, t and x, 
FiV'it) > x,J\fi{ex,t) =0} 



\ sup Sin> x{l-e{ri-di)),Mi{ex,t) = o\. 

^n<N^(t) J 



',Nf{t) 
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Proof. We have 

F{V''{t)>x,Miiex,t)=0} 

<F{vf^{t)>x,Miiex,t) = 0} 



<p|(ri-(i,)(BiA[o + (l-B,)Ai 
+ sup Sin> x,J\fi{ex,t) = 0> 

n<N^(t) i 

<f\ sup Sin>x{l-e{ri-di)),Afi{ex,t) = o\. 



i<N^{t) 

The last inequahty follows from the fact that A[q and Ajo must be smaller 
than ex if A/"j(ex,t) = 0. D 

To obtain upper bounds for probabilities as in Lemma 6.1, we will fre- 
quently apply the following key lemma, which is a trivial modification of 
Lemma 3 in [34]. 

Lemma 6.2. Let S„ = Xi + • • • + X„ be a random walk with i.i.d. step 
sizes such that E{Xi} < and E{(X^)p} < oo for some p> 1. Then, for 
any /3 < oo, there exists an e* > and a function </>(•) G 7^-/3 such that, for 

ee(0,e*], 

P{S„ > X I Xj < ex,j = 1, . . . ,n} < 0(x), 
for all n and all x. 

Note that if Xj can be represented as the difference of two nonnegative 
independent random variables X^ and X^, then the lemma remains valid if 
the Xj's are replaced by Xj. 

The final preparatory lemma is a simple consequence of Corollary 4.2, 
which will be used several times in combination with Lemma 6.2 that prob- 
abilities of certain events are of o(P{V'^ > x}). Define 

N N 

P{x) := n 1P{A5 > x} G 7^_^, ^^ := ^(z., - 1). 



Lemma 6.3. 



limsup T—-— ^ < oo. 

^^oo^ P{Vc > x} 
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We now show that, with overwhelming probabihty (as x -^ oo), the rare 
event {V^ > x} occurs as foUows: 

(i) The process {A{0,t) — ct} reaches level x before time Mx for some 
large M. 

(ii) Up to time Mx, each flow generates exactly one long On period, that 
is, Mi{ex,t) = 1 iov i = l,...,N. 



Proposition 6.1. 



^. ^. . ^F{V%Mx) > x} ^ 
lim hmmt = 1. 

Af^oo x->oo P{V^ > x} 



Proof. By definition, 
FiV >x}= p(sup{A(0, t) -ct}>x 



t>0 



<p( sup {A{0,t) - ct} >x\+f\ sup {A{0,t)-ct}>x 

[o<t<Mx } [t>Mx 

= FiViMx) >x} + f\ sup {A{0, t)-ct}>x\. 

[t>Mx J 

Thus, it suffices to show 

,. ,. F{supt>MAMO,t) -ct}>x} 
lim lim sup ~ ^,, ^ ^ = (J, 

which, however, follows directly from Lemma 4.5. D 

Now suppose that the workload reaches level x. By the previous propo- 
sition, we may assume that this occurs before time Mx (for M sufficiently 
large). The next two propositions show that we may restrict the attention 
to a scenario where each flow initiates exactly one long On period before 
time Mx. 

The first proposition indicates that each flow has at least one long On 
period. 

Proposition 6.2. For all i and all M, there exists an e* > such that, 
for all £ £ (0,e*], 

F{V^{Mx) > x,J\fi{ex, Mx) = 0} = o(P{V^ > x}), 

as X ^>- oo. 
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Proof. Define A^f (t) := max{n : ^2]=! ^ij <*} + !• Note that N^^{t) < 

Nf{t). 

Using Lemma 6.1, taking t = Mx, 

F{Y^{Mx) > x,Mi{ex, Mx) = 0} 



<p( sup Sin>xil-e{ri-di)),Mi{ex,Mx)=o] 

'■-n<Nf{Mx) J 

<p| sup Sin>x{l-e{ri-di))\Mi{ex,Mx)=o\ 

= F< sup Sin> x{l- e{ri-di))\Aij <ex, 

j = l,...,Nf{Mx) 
= W\ sup Sin>x{l-e{ri-di))\Aij<ex,j>l\ 

<pi sup Sin>x{l-£{ri-di))\Aij<ex,j>l\ 

U<iVf(Mx) J 

= ¥< sup Sin>x{l-£{ri-di))\Aij <£x, 

U<iVf (Mx) 

j = l,...,N^{Mx) 

<f\ sup Sin>x{l-£{ri-di))\Aij<£X,j>l\ 
in<M2X J 

+ F{N^{Mx)>M2x}. 

The second term decays exponentially fast in x if M2 > XiM. The first term 
can be bounded by 

M2X 

Y^ F{Sim>x{l-£{ri-di))\Aij <ex,j = l,...,m}. 

m=l 

According to Lemma 6.2, there exists an e* > and a function (j)(-) G 7^_/3 
with (3 > fi + l, such that, for e € (0, e*] , the last quantity is upper bounded 
by M2xcj){x). The latter function is regularly varying of index 1 — /3 < — /U. 
Invoking Lemma 6.3 then completes the proof. D 

The next proposition shows that each flow has at most one long On period. 
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Proposition 6.3. For all i, all M and all e > 0, 

¥{V^{Mx) > x,Mi{ex,Mx) > 2} = o(P{V" > x}), 
as X ^ oo . 

Proof. Without loss of generality we may take i = 1. By Proposition 6.2 
it suffices to show that 

F{V%Mx) > x,Mi{ex,Mx) > 2,J\fi{ex,Mx) > l,i > 2} = o(P{V^ > x}). 

Note that the left-hand side is bounded by 

N 

¥{Afi{ex,Mx)>2}Y[F{J\fi{ex,Mx) > 1}. 

Thus, invoking Lemma 6.3, it suffices to show that: 

(i) ¥{J\fi{ex, Mx) > 1} is bounded by a function which is regularly vary- 
ing of index 1 — Ui. 

(ii) ¥{Mi{ex,Mx) > 2} = o(F{Mi{ex,Mx) > 1}). 

We will prove both assertions for i = 1. For assertion (i), note that 

P{7Vi(ex,Mx)>l} 

< piF{Al > ex} + P{#{j G {1, . . . , N^^iMx)} : Ay > ex} > 1}. 

The first term is in TZi-,/^ . By conditioning upon N^ {Mx), the second term 
can be bounded by ¥j{Ni {Mx)}¥{A.i > ex}, which is also regularly varying 
of index 1 — i^i . To prove assertion (ii) , note that 

¥{Mi{ex,Mx) >2} 

<Pi¥{Al>ex}F{Mii{ex,oo),{0,Mx]) > 1} 

+ F{Mi{{ex, oo), (0, Mx]) > 2}. 

Using F{Afi{{ex,oo),{0,Mx]) > 1} < F{Afi{ex,Mx) > 1} and assertion (i), 
it follows that the first term is of o{F{Afi{ex,Mx) > 1}). To bound the 
second term, condition (again) on N\^{Mx). This yields 

P{A/'i((ex,oo), (0,Mx)) > 2} <E{iVf (Mx)^}P{Ai > ex}2. 

Finally, note that E{iV|^(Mx)^} is quadratic in x for x — > oo. D 

We have now shown that, with overwhelming probability, each flow con- 
tributes to a large value of sup^>o{j4(0,t) — ct} through exactly one long 
On period. We thus proceed with the second statement (as indicated at the 
beginning of this section), implying that apart from these long On periods, 
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the flows show typical behavior. In order to formalize that statement, we 
need to introduce some notation. Define 

T{y) := inf{i > : ^(0, t)-ct = y} 

as the first time at which the process {A{0,t) — ct} reaches level y. 

For fixed e > and x, let Ts^i{£x) and Tf^i{ex) be the respective starting 
and finishing times of the first On period of flow i exceeding length ex. 
Denote 

Ts(ex) := max Tsi(ex) 
^ ^ i=h...,N '^ ' 



and 



Tf{ex) := .jnin^T/,i(ex). 



Note that all flows are in the middle of their long On periods between 
times Ts{£x) and Tf{ex). We will show that the fluctuations of the process 
{A{0,t) — ct} away from the mean before time Ts{£x) and after time Tf{£x) 
can be neglected. 

More formally, the next two propositions show that, given that the work- 
load reaches level x before time Mx, there exists for any small 5 > an e^ 
such that, for all e G (0,65), 

Tsiex) < t{6x) < r((l — 6)x) < Tf{ex). 

Thus, the workload remains small up to time Ts{ex), and reaches a level 
close to X before time Tf{ex), as depicted in Figure 2. 

The first proposition indicates that it is most unlikely that the process 
{A{0,t) — ct} reaches level 6x before time Ts{ex). 

Proposition 6.4. For any 6 > 0, there exists an e* > such that, for 
all e e (0,e*], 

F{t{6x) < Ts{ex)} = o(P{V^ > x}). 
Proof. For compactness, denote Ts = Ts{ex), Ts^i = Ts^i{ex). Then 

N 

F{Ti5x) < Ts} = ny"irs) > Sx} <Y.^{V''iTs,^) > 6x}. 

i=\ 

We bound each term in the last summation. 

Define Ni{sx) := N^{t~j) as the number of On periods initiated by fiow i 
before the first On period exceeding length ex. Note that Ni{ex) + 1 is 
geometrically distributed with parameter P{Aj > ex}. 
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T((l- 5)x) T^(ex) ^ 



Fig. 2. Typical path to overflow. 

Using Lemma 6.1, taking t = Tg^i, 
nV"{rs,^) > 6x} 

= ny%rs,i) > 6x,J\fii{ex, oo), [0,r,,,)) = 0} 

<P<^ sup Sn>x{5-e{ri-di)),Aij<£x,j = l,...,Ni{£x)\ 

in<Ni{ex) J 

oo 

m=l 

X P<^ sup S„ >x{6 -e{ri -di)),Aij <£x,j = l,...,m[ 

oo 



m=l 



X P<^ sup Sn > x{6 -e{ri - di)) \Aij <ex,j = l,...,mL 

According to Lemma 6.2, there exists an e* > and a function (j){-) € TZ-/3 
with /? > 2i/ + 1, such that, for e G (0, e*), the last quantity is upper bounded 

by 

f){x)F{Ai < ex} 



E{Ni{ex)}(t>{x) 



P{ Ai > ex} ' 
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which is regularly varying of index i/j — /3</i + l — {2fj, + 1) = — /i. 
Invoking Lemma 6.3 then completes the proof. □ 

The next proposition shows that, given that the process {A{Q,t) — ct} 
reaches level x before time Mx, most probably level (1 — 5)x is crossed 
before time Tf{ex). 

Proposition 6.5. For any 6 > and M < oo, there exists an e* > 
such that, for all e G (0,e*), 

P{r((l - 6)x) > Tf{ex),V''{Mx) > x} = o(P{V^ > x}). 

Proof. For conciseness, denote tj = Tf{ex), Tf^i = Tf^i{£x). By Propo- 
sitions 6.2 and 6.3, it suffices to show that 

P{r((l - 6)x) > Tf,V''{Mx) > x,J\fi{ex, Mx) = 1 for alH = 1, . . . , iV} 

= o(P{V">x}). 

Note that 

P{r((l - 5)x) > Tf, ViMx) > x,Mi{ex,Mx) = 1 for alH = 1, . . . , iV} 

= P{V^'(t/) > (1 - 6)x, V'iMx) > X, 

Mi{ex,Mx) = l for a\li = l,...,N} 

N 

< ^P{y"(T/,0 > (1 - 5)x, V^{Mx) > x,Ni{ex,Mx) = 1}. 

As before, we bound each term in the last summation: 
P{y^(r/,i) > (1 - 5)x, V^{Mx) > x,Mi{ex, Mx) = 1} 

<p| sup {A{0,t)-ct}<{l-6)x, sup {A{0,t)-ct}>x, 

lo<t<Tf^i 0<t<Mx 

Mi{iex,(X)),{Tf^i,Mx]) = 

<P( sup {A{Tf^i,t)-cit-Tf^i)}>6x, 
KTf^i<t<Mx 

7Vi((ex,oo),(T/,i,Mx]) = 
<p| sup {Ai{Tf^^,t)-di{t-Tf^i)]>5x, 

yTf^i<t<Mx 

A/'i((ex,oo),(T/,i,Mx]) = 
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The first inequality follows from the definitions. The second inequality fol- 
lows from properties of the sup operator, and the last inequality is obtained 
by assuming that all flows but i are On between times Tf^i and Mx. 
Note that the last probability is upper bounded by 

P\ sup Sn-SNi(ex)+l>^X, 

^Ni{sx)+2<n<N^{Mx) 

Aj <ex,Ni (ex) +2<j< Nf{Mx) \ . 

The latter probability can be upper bounded by a function which is regularly 
varying of index — /3 < —fi in a similar fashion as in the proof of Propositions 
6.2 and 6.4. 

The proof is complete by invoking Lemma 6.3. D 

Propositions 6.4 and 6.5 may be used to obtain the following result. 

Corollary 6.1. If Aj{-) eTZ for all j = 1,...,N, ^/len P{V^ > x} G 

in. 

The above result suffices to prove the reduced-load equivalence. (See Sec- 
tion 5, in particular Proposition 5.1, for the details.) However, determining 
the exact asymptotic behavior of PjV'^ > x} requires further analysis, to 
be found in Sections 6.3 and 6.4. In particular, the analysis in Section 6.4 
will lead to a sharper version of Corollary 6.1, showing that PjV^ > x} (^TZ 
(which is a strict subset of TTZ) . 

Nevertheless, we sketch a direct proof of Corollary 6.1 which we believe 
is of independent interest. For the formal proof details we refer to [41]. 

Sketch of proof of Corollary 6.1. The idea of the proof is as 
follows. If V^ > X, then Propositions 6.4 and 6.5 show that the process 
{A{0,t) — ct} reaches the level (1 — 5)x after all flows have been On for at 
least ^ ~_ , ^ time units. Since Aj{-) G 7?. C ITZ for all j = 1, . . . , A^, with high 
probability, all flows remain On for at least ^r^ more time units. This yields 

limliminf PIV^ >(l + 6)x I V^ > x| = 1, 

implying the desired statement (by definition). D 

6.3. Proof of Theorem 6.1. In this section we give a proof of Theorem 6. 1 . 
First we consolidate the key results from the previous section in the following 
theorem. 
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Theorem 6.2. For any 6 > 0, there exists an e* > such that, for all 
eG(0,e*), 

F{A{0,Tf{ex)) - cTf{ex) > x} < PfV^ > x} 

< F{A{0, Tfiex)) - cTfiex) > (1 - 5)x}. 

Proof. The lower bound is trivial. The upper bound follows from Propo- 
sitions 6.1, 6.4 and 6.5. D 

In order to obtain tight bounds for the probabilities in Theorem 6.2, we 
condition upon Tg i for all i. Hence, for any Jq C J ^ define the event Djg{£x) 
by 

Djg{ex) := {Ts^i{ex) = for all i G Jo;Ts^i{ex) > for all i ^ Jo}. 

The event Dj^ (ex) implies that the flows i ^ Jo started their long On period 
before time (remember that we consider the system in stationarity). The 
flows i G J7i start their long On period at a later time epoch. 

Denote Pj^oj-} = P{- | Djq{£x)}. The following two lemmas will be useful 
for providing tight upper and lower bounds for the probabilities in Theo- 
rem 6.2. 

Lemma 6.4 (Upper bound). For any 5 > 0, there exists an es > such 
that, for all e G (0,6^), 

P^„{^(0,r/(ex)) - cTfiex) > (1 - 5)x} [] HK > ^x} 
<Pjoi{l-S)x)l[pi, 

with Pjq{{1 — 6)x) as in (6.1). 

Lemma 6.5 (Lower bound). There exists an e > such that 

Fj,{A{0,Tf{ex)) - CTfiex) > x} J] P{A[ > ex} 

ieJo 



r^Pjoix) n Pi' 



with Pjg{x) as in (6.1). 

The proofs of these lemmas are quite technical, and are deferred to Ap- 
pendices A and B. A brief sketch of the proofs is given at the end of this 
section. 

We now have gathered all the ingredients for the proof of Theorem 6.1. 



FLUID QUEUES WITH HEAVY TAILS 37 

Proof of Theorem 6.1. The lower bound in Theorem 6.2 may be 
written as 

F{A{0,Tf{ex)) - cTf{ex) > x} 

= Yl ^Jo{MO,Tf{ex)) - cTfiex) > x}F{Dj,{ex)}. 

JoC{l,...,N} 

Note that 

F{Dj,{ex)}r^ Ylp^F{Al>ex}. 

i&Jo 

Using Lemma 6.5, we then obtain 

F{A{0,Tf{ex)) - CTfiex) > x} 



Vi=l / JoC{l,...,N} 



Similarly, using Lemma 6.4, 

Fj,{A{0,Tfiex)) - CTfiex) > (1 - 6)x} 

V=l / JoC{l,...,N} 
Theorem 6.2 then gives 



fn^.) E pjoi^) 

Vj=l / JoC{l,...,N} 

<F{V'>x} 

N 



\j=l I JoC{l,...,Af} 



which implies Theorem 6.1, since Pj^ix) G 7^ as will be shown in Theo- 
rem 6.3. 

D 

In preparation for the proofs of Lemmas 6.4 and 6.5, we give a convenient 
representation for yl(0,rj) — cTf under the event Dj^^iex). 



as 



Lemma 6.6. Under the event Dj^iex), A(0, Tf) — cTf can be represented 

AiO,Tf)-CTf 



min < min F j , min G j 
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where J\ = J\Jq. The random variables F j and Gj are given by 
Fi = ir-c)Aliex) 



keJi 



BkAl{ex) + {l-Bk)[Ak{ex)+\Jl]+ ^ Vkj 



Gj = {r-c)Ai{ex 
+ ir-c) 



N,{ex) 



di 



Niiex) 

(1-B,)U[+ Y. U., 





- 


Nk{ex) 


>; n 


(1- 


-Bk)vi+ >; Uk, 


keJi\{i} 


- 


J=l 



Here Ai{ex) = Aj|Aj > ex, A[(ex) = A[|A[ > ex, Aij{ex) = Aij\Aij < ex 
and A[(ex) = A[|A[ < ex. 

Proof. Under the event Dj^^ex), the random variables Tg^i, i € J'l, can 
be represented as 



Niiex) 

Ts,i = BiAl{ex) + {l-Bi)[Ul + Ai{ex)]+ ^ [U^^- + A^^- (ex) 



ieJi- 



Combined with the identities 
Ai{0,Ts,i) =ri 



Niiex) 



BiA'i{ex) + {l-Bi)Ai{ex)+ J2 Aij(ex) 



Tf = min< min AJ'(ex),min{Ai(ex) +Tsi} 

the representation for A{0,Tf) — ctj then easily follows. D 

We now give a brief sketch of the proofs of Lemmas 6.4 and 6.5. Both rely 
on the above representation for j4(0,tj) — cTf in terms of the variables Fj 
and Gj. The proofs of the lemmas have a similar structure. 

(i) The expressions for Fj and Gj are quite complicated, so an attempt 
to obtain the exact joint distribution does not seem promising. Therefore, 
the first step is to show that all random variables Ajj(ex) and Ujj can be 
replaced by their means. 
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(ii) The above point indicates that Fj and Gj may be approximated as 
follows: 

Fi^{r- c)Al{ex) + ^ nK{Vk}Nk{ex), 
keJi 

Gi « (r - c)A,(ex) + [(r - c)E{AJ - d,E{^i}]Ni{ex) 

- Y. nE{lJk}Nk{ex). 
keJi\{i} 

It will be useful to keep these approximations in mind. The formulas in 
Appendixes A and B look much more cumbersome by the appearance of 
many additional but small constants. 

(iii) The only random variables appearing in the above expressions are 
Ai{£x), B[(ex) and Aj(ex), of which the distributions are known. What 
thus remains is a straightforward computation. 

The first point causes the most technical difficulties. It requires a separate 
treatment in the proofs of Lemmas 6.4 and 6.5. Details may be found in the 
Appendices. 

6.4. Asymptotic behavior of Pjg{x) and ¥{V'^ > x}. In this section we 
give an asymptotic characterization of Pjg{x), which may be useful for fur- 
ther analysis. In particular, we establish that Pjq{x) and PjV^ > x} are 
both regularly varying. 

Assume that Jq is a proper subset of J , observing 



Pj{x) = X{Aa', 



X 

> — 



r 



For every set J^o, define the |^i|- vector g by 

i - Pj 



9 , 

\r-c J j^j^ 

Let G be a (square) matrix with identical rows <?, and let G .= G — I, with 
/ the identity matrix of dimension \J'i\- 

It can easily be shown that G is invertible; denote its inverse by H. 
A straightforward computation yields H = }_i G — /, with e = (l,...,l), 

which implies that gH = ^_^ g- A further straightforward computation shows 

1^1=65-1. 

Define y = {yi)j^ and dy = Iliej-i '^yi-- Then we may write 

X I ^^\{A^.>{Gy\ + ^\\{AK>9y + ^\dy. 
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If we integrate w.r.t. z := Gy (note that (5 is a positive matrix), then we 
obtain [defining Aj^ = (A[)iej7j 

1 -^Ar>^,^Gj; 



Al--^>-^g(A^j^-e-^_],keJ,\. 



1 i 7 

I r — c 

X 1 / X 

> 9[ At, - e 

r — c eg — 1 \ r — Cy 

We conclude that Pjg (x) can be written (up to a constant) as the probabihty 
that (Apjgj- belongs to a certain set. We now show that Pj-q(x) is regularly 
varying of index —fj, [recall that fj, = J2i=i{'^i ~ !)]■ If -Aj is regularly varying 
of index —Vi < —1, then it is well known and elementary to show that 



as X ^ oo. Let Zj be a random variable with the above limiting distribu- 
tion, with 7 = ^— such that the Zj, i G J7i, are independent. The above 
computations are summarized in the following theorem. 

Theorem 6.3. 

N 



with Kj = 1 and 



i=l 

- 1 I eg -I J 



65 -1 

if Jo is a proper subset of J . In particular^ Pjoi^) ^^ regularly varying of 
index —fi. 

Combining Theorems 6.1 and 6.3, we obtain: 

Theorem 6.4. 

N 



( ^ 

{N^>x}^>^^pAa\>^V 
1=1 
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with 

fi= J2 ^Jo- 
JoQ{l,...,N} 

In particular, PjV^ > x} is regularly varying of index —fi. 

The above theorem is used in proving the reduced-load equivalence (see 
Section 5), and may be potentially useful for computational purposes. 

In particular, in the case of two On-Off flows, the computation of k is as 
difficult as the computation of ki and K2- Using the probabilistic interpre- 
tation of these constants readily leads to an integral expression, which can 
be solved explicitly when both vi and 1^2 are integer-valued. We omit the 
details. 

7. K heterogeneous classes: proofs. In this section we provide the proofs 
of the results in Section 3.5 for the case with K heterogeneous classes of On- 
Off flows. In particular, we present a proof of Theorem 3.3. 

We start with the regime where we first let x ^ 00 and then n ^ 00. For 
every n we have, using Theorem 3.2, 

P{V(")>nx} .„^ 

hm — !^ — i = -/iW^ 

x^oD log X 

with /i^"^' denoting the optimal value of the criterion function of the associ- 
ated knapsack problem. It thus remains to be shown that 

(7.1) lim ^— = 1. 

First observe that the optimal value of the continuous relaxation of the 
knapsack problem is n/i, yielding a lower bound for n^"^'. On the other hand, 
the continuous relaxation may be used to construct a feasible solution of the 
knapsack problem. Take (use the notation of Section 3.5) qk = n^ = na^ for 
k <l, qk = Uk = for k> I, and qi = |n/| -|- 1. This is a feasible solution with 
a value at most nfi + ^i, giving an upper bound for fi'^'^'. In conclusion, we 
have 

nfi < fi^"'' <nfi + 'ji, 

from which (7.1) directly follows. 

We now turn to the regime where we first let n ^ 00 and then x — > 00 
(i.e., the many-sources regime). Define the "decay rate" 

I(x) := - lim ilogP{V(") >nx}. 

It needs to be shown that /(x) ~ /ilogx as x — s- 00. 
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The above decay rate equals (8, page 300) 

I{x) = mi sup 1 9 {x + t)-f2 afclogIE{e^^^W} j , 

with Ak{t) := Ak{0,t) representing the amount of traffic generated by a 
single class-A; flow in a time interval of length t in steady state. Replacing 6 
by 0{logt)/t, we obtain an alternative variational problem, 

inf logi- Jt(- + 1 

(7.2) 



where Jt{x) := sup 6x — >^ a^ 



/ #, logE{e^(i°g*)^fcW/n\ 

for x € {0,J2k=i^kfk)- The latter variational problem allows direct asymp- 
totic analysis (x — >■ oo) as in [29], which yields Theorem 7.1. 

First, however, we state an auxiliary lemma. Recall that a^ = J2m=i '^mfm + 
J2m=k^mPm, and that the K classes are indexed in nondecreasing order of 
the ratios jk = (i^fc - l)/(?^fc - Pk)- 

Lemma 7.1. For9>0, 

logE{e^('°g*)^.WA} 

hm — = max|6'pfc, 6'rfc - i/^ + 1}, 

t-^oo logr 

so that the cumulant function of the superposition is piecewise linear: 
V Ofc lim — = y] «fc max{6'pfc, 9rk -i^k + !}• 

Further, 

l{x)-l K 

C^-^) t^"^^^^^ " ^'(^)^ ~ ^ ^k{ll{x)rk - ffe + 1) - Yl 0,kll{x)Pk, 

^°° fc=l fc=i(a:) 

/or X G (0,X]fc=i '^fc'^fc)' where l{x) is such that x G ('7/(x)-i)<7i(a:))- 
The function linit^oo Jti') is increasing. 

The proof of the above lemma is analogous to that of Theorem 3.6 and 
Lemma 3.7 of 



Theorem 7.1 (Large-buffer asymptotics) . 

hm = fi, 

x^oa log X 

with /U = Ei=i«fc(i^fc - 1) + (1 - cri-i)^i and I :=?(!)• 
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Proof. The proof consists of deriving an upper bound and a lower 
bound which asymptotically coincide. 

(Upper bound.) Using the representation (7.2), 

r -^(^) r ■ f ^°S* J fx \ 
limsup- = limsupmi JJ — I" 1 • 

x^oo logX x^oo OOlogX \t J 

Substituting t = x/s, s G {0,J2k=i o-k^k — 1); to obtain an upper bound, and 
using (7.3), 

lim sup inf Jt { — hi 

x^oo i>OlogX V* 

< limsup f ' V ,./, g + 1 

x^oo log X ' 

log(a;/s) f , -,\ 

< limsup — limsup J^/j,(s + Ij 

x— »oo log X X— >oo 

< limsup J2./s(s + 1) 

X— »oo 

= 7/{s+i)(s + l) 

lis+l)-l K 

- 51 {<^kll{s+l)rk-T^k + '^)- XI (^kll{s+l)Pk- 
k=l k=l{s+l) 

The above inequality holds for any s G (0,X]fe=i'^fe''A: — !)• According to 
Lemma 7.1, the last term is increasing in s + 1. Letting s J, to obtain the 
sharpest possible upper bound, we obtain 

I(x) '"^ ^ 

limsup <7i-y] akijiTk - I'fc + 1) - V a^'jiPk = P- 

™ logx ^^ ^^ 

(Lower bound.) Using the representation (7.2), and taking 6 = ^i, we 
obtain the lower bound 



I{^) = inf logt • sup (9 ( - + 1 ) - X Ofc 



*>o e \ \t J t^i log* 



K 



logE{e^('°s*)^fcW/*}~ 



>inflogt. 7.7 + 1 HE 



x 



K 



log E{e'>'' ('°s i)Ak (t)/t 1 \ 

The optimizing value of t in the above variational problem is at least linear 
in x, for large x. Formally, there exists a d such that the above infimum need 
be taken only over t > dx, for large x. This may be proven analogously to 
case (iii) of [14], page 258. Thus, 

I{x) > m£logt . (^7. [- + !)- E «. ^^ j . 
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Using (7.3), we find that for any e > 0, and x large enough, we have, for all 
t > dx, 

^ log E{e^'i°s*^^- (*)/*} ^^^ ,^ 

Vflfc ] — : < (l + e) > afemax{7//9fc,7zrfc-ffc + l}. 



Thus, 



lim inf 

x^oo logX 



> lim inf inf :; 7/1 — hi 

x^oo t>dx log x\ \t 



K \ 

(1 + e) ^ Ofc max{7ipfc,7zrfc -Uk + l}\ 
k=l J 



> lim inf inf inf 7M — hi 

x^oo t>dx log X t>dx \ \t 



(1 + e) ^ ttfc max{7;pfc, 7irfc - i/fc + 1} | 
fc=i / 



ii" 



> 7i - (1 +e) ^ aA;max{7;pfc,7irfc -i^k + 1}- 



fc=i 
Letting e j 0, we obtain 

lim inf : 



x^co log X 

K 

>7i-^ak max{7//3fc, 7;rfc - i^^ + 1} 
fe=i 

K 

= 11-^ a-kiliPk + max{0, ^i{rk - pk) -i^k + 1}) 
k=l 

K 

= 11 -Yl o-killPk + max{0, {-fi - -fk){rk - Pk)}) 
fc=i 

= 7z - X! o.k{lirk - i^fc + 1) - X! OfcTi/Ofc 
fc=i fc=i 

= p. U 
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As shown above, Theorem 3.3 impHes that the Umits x ^ oo and n ^ oo 
commute, as long as one considers "rough" (i.e., logarithmic) asymptotics. 
However, in case of "more refined" asymptotics, the limits do not necessarily 
commute. This may be seen as follows. Consider the case of n homogeneous 
On-Off flows with Pareto(z/) distributed On periods. In [28], it is proven 
that 



lim lim -logPlV'-"'' > nxj + (z^ - 1) ^ logfxlogj;) = i/, 

a;— >oon— >oo fi \y — n I 

for some constant H £ (0, oo). Now reverse the limits. Denote by A;„ the 
number of flows sending at peak rate in the reduced-load approximation (in 
the notation of Section 3.4, we have kn = N*): 



k 



nc — np 

n 

r — p 



Now with Theorem 3.1, we have, for any flnite n and a; — > oo, 

P{V(") > nx} ~ /(njx^^'^-i)^", 
for some function /(•). Hence, 

lim - log P{ V(") > nx} + {v-l){ - — - ) log(j; log x) 
x^oo n \r — pj 

= log/(n) + lim (z^ — 1) — I logx — (z/ — 1) log log x. 

x-»oo \ n r — pj r — p 

Since this limit does not exist in M, we conclude that the limits do not 
necessarily commute. 

8. Concluding remarks. We have characterized the asymptotic behavior 
of the workload distribution in a fluid queue fed by multiple heavy-tailed 
On-Off flows. The results extend previous work, like the bounds derived 
in [15], and the exact asymptotics in [9, 22] which rely on strong peak-rate 
conditions. As a by-product, the proofs lead to several important insights 
like the extension of the reduced-load equivalence established in [1] (see 
Section 5), and a detailed understanding of the typical overflow behavior 
(see Section 6). In the analysis, we excluded the case where the drift may 
be zero during the path to overflow (see Section 3.1 for a brief discussion), 
which appears particularly interesting from a theoretical perspective. 

There are several other interesting topics for further research. We expect 
that the methodology of Section 6 is also suitable to study other similar 
problems, such as fluid queues with M/G/oo input, multiserver queues, and 
Generalized Processor Sharing queues. A further avenue for research is the 
extension of the results to the case of On-Off flows with more general subex- 
ponential On periods, for example, Weibull. Partial results in [1] indicate 
that the typical overflow behavior may then actually be quite different. 



46 B. ZWART, S. BORST AND M. MANDJES 

APPENDIX A 
Proof of Lemma 6.4. 

Lemma 6.4 (Upper bound). For any 6 > 0, there exists an es > such 
that, for all e G (0,65), 

P^„{^(0,rj(ex)) - cTfiex) > (1 - 5)x} \{ P{A[ > ex} 
<Pj,{{l-6)x)\{p,, 

with PJ^^{{\ — 6)x) as in (6.1). 

Proof. As mentioned earlier, the first step is to replace all random 
variables Ajj and Ujj by their means. Let 6 and 5 be two |j7i| -vectors, of 
which the elements are positive but arbitrarily small. Note that, for fixed 
Jo, 

Fi<{r- c)Al{ex) - Y, nNk{ex)[E{Vk} - 4] 

k&Ji 
Nkiex) 

+ ^rfc ^ [E{Ufc}-4-Ufc,], 
keJi j=i 

Gi < (r - c)Ai{ex) + (r - c)ex + (r - c)Ni{ex)[E{Ai} + 6i] 

N,{ex) 

+ {r-c) Y. [Aij{ex)-E{A,}-~5,]-diN,{ex)[K{lJ,}-~6i] 
j=i 

Niiex) 

+ d, J2 m^^}-Si-Vij] J2 rkNk{ex)[E{\Jk}-6k] 
i=i keji\{i} 

Nkiex) 

keJi\{i} j=i 
Define the event £'1(7,5, 5, e,x) by 

I J2 [E{U i} -mm{6ji} -Vij] <jx/{2r),i e jA 

r Niiex) 

ul Y [A^j{ex)-E{A^}-mm{5i,6^}] 
[ j=i 
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< ^x/{2r) — (r — c)ex,i £ Ji>. 

A straightforward application of Lemma 6.2 (analogously to the proofs of 
Propositions 6.2, 6.4 and 6.5) shows that, for any -^,5,6 > 0, there exists an 
e* > such that, for all e G (0,e*], 

(A.l) ¥j„{Ei{j,6,~6,e,xr} = o{P{x)), 

as X ^ oo with P{x) = nj=i]P{A^ > x}, as defined earlier. 

From (A.l) and Lemma 6.6, we conclude that, using the upper bounds 
for Fj and Gj, 

¥j,{A{0,Tf)-cTf>{l-6)x} 

= Fj,{A{0,Tf) - CTf > {I - 6)x;E^{-,,lle,xy} 

+ Fj,{A{0,Tf) - cTf > {I - 6)x;E^{^,lS,e,x)} 

<p|(r-c)A[(ex) 

- E rkNk{ex)[E{Uk} - 4] > (1 - 7 " S)x,i& Jo; 
fceJi 

(r - c)Mex) + (r - c)iV, (ex) [E{AJ + 5^] 

-diNi{ex)[E{V^}-~6i] 

- Yl rkNkiex)[E{Uk}-5k]>il--/-6)x,ieJi\ 
keJi\{i} ) 

+ o{P{x)). 
The last probability equals [condition on Ni{ex), i € J'l] 

ni>l,i£Ji \ieJi ) 

xpJ (r-c)A[(ex) 

- Y^ rfc[E{Ufc}-4]nfc>(l-7-'^)a;,iG Jo; 
(r - c)Ai(ex) + (r - c)[E{A,} + l^Ui 
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^ rfc[E{Ufe} - ~5k\nk > (1 - 7 " ^)x.i^ j\ 
J^\H^ ) 



k&Ji\{i\ 

Deconditioning upon A^ and A[ (i.e., dividing by IliG Jo ^i-'^i -^ ^^i Ilig Ji IPjAj > 
ex}), and noting that P{A'j(ex) = rii} < P{Aj > ex}, we obtain that 

¥j,{AiO,Tf)-cTf>{l-5)x} 

X Yl P{A[ > ex} 
ieJo 

is upper bounded by [up to o{P{x))] 



ni>0,ieJi \ieJo ^ 

+ J2 r-fe[E{Ufc}-4]nfc 



k£ji 



Yl P<(r-c)Ai>(l-7-5)x 
ieJi I 



X 

ieJi 



+ [diK{Ui} - (r - c)E{Ai} - ri6i]ni 



+ Yl rk[E{lJk}-~6k]nX 
ke.7AU\ ) 



fceJi\{i} 

It is important to note that this expression is independent of e. 

Since ah probabihties appearing in the right-hand side are decreasing 
functions of Ui (for 6 and 5 small enough), the latter term is bounded by 
[with y := {yi)i(zji and dy := UieJi ^V'^ 

[ p| (r - c)A[ >{l-^-5)x+J2 ^fc[E{Ufc} - h]yk \ 
Y[ P< (r - c)A, > (1 - 7 - (5)x 

iGJi I 



(A.2) 

+ [diE{U,} - (r - c)E{A,} - nSijyi 



+ Y nm^k}-Sk]yk]dy. 

keJAU} ) 



k&Ji\{i} 

We will rewrite this expression in terms of Pjf^[x). Apply the change of 
variables Zi := ?/i/(E{Ai} + E{U,J). Redefine 5i := (5i/(E{AJ + E{Ui}) and 
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similarly ~5i := 6i/{E{Ai} + EjUj). Note that 



E{Aa + E{U.J E{Ai} 'E{Ai} + E{Ui} 

Then we obtain that (A. 2) equals 



n 



Pi 



..e^.JE{AJ^ 



/ J]pj(r-c)A[>(l-7-5):E+ ^(r.-p,-^)^^ 
X n IP"] (^ - c) Ai > (1 - 7 - '^)2; + (di - /Oi - 5i)2i 

+ H {rk- Pk- h)zk \ dz. 



n 



1 l-l-Q- QT^ff in + o rr V a "I" o wt v "I" 'y ■ —1 i-1 1 

fi pi 

ri - Pi Pi 



If we take 5,- = ' ^' J,- and integrate w.r.t. z,- ''' ^' — ^, then we obtain 



^iGJl''*"^*"'^^^^^*^' 



X / []pj(r--c)A[>(l-7-5)x+ ^(rfc-pfc)^4 
X n 'P'l (^ - c)Ai > (1 - 7 - 5)x 

+ {di - pi)zi + Y^ (rfc - pk)zk > dz 
keJi\{i} J 

ieJi ^^ P' ^* 

Together with the fact that Pjq{-) is regularly varying (see below), this 
completes the proof of the upper bound after dividing by Pjq{x), letting 
X -^ (X), and noting that 6, 6 and 7 may be chosen arbitrarily small. D 

APPENDIX B 

Lemma 6.5 (Lower bound). There exists an e > such that 
Fj,{A{0, Tfiex)) - cTfiex) > x} H P{A[ > ex} > Pj,{x) J] Pu 

i&Jo i&Jl 
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with Pjq{x) as in (6.1). 

Proof. As in Appendix A, the first step is to replace the random vari- 
ables Aj(ex) and Uj by their means. Adding and subtracting appropriate 
means, it is easy to see that, for fixed j7o, 

Fi = (r - c)A\{ex) - ^ rkNk{ex)[¥.{\]k] + 4] 

+ ^ rfc ^ [E{Ufc}-Ufcj+4] 

- Y. rfc[BfcA^(ex) + (l-Bfc)(Afc(ex)+UD], 

k&Ji 

G, = (r - c)ki{ex) + (r - c)[BiA[(ex) + (1 - B,)A,(ex)] 
-d,(l-B,OU[- Y. ^fc(l-Bfc)Ul 

k(^Ji\{i} 
N,(ex) 

+ {r-c) Y [Aiji£x)-E{A,} + ~6,] 
i=i 

+ ir-c)Ni{ex)[E{A^}-5i] 

-d^N^{ex)[E{Ui} + ~6i] 

N,{ex) 

+ d^ Y []E{U,}-U,,+5i] 

- ^ rfciVfc(erc)[E{Ufc} + 4] 
keJi\{i} 

Nkiex) 

+ 5^ rfc ^ [E{Ufe}-U,, + 4]. 

fc6j"l\{i} J = l 

Define the event £'2(7,^, (5,e,x) by 

j Y [E{Ui}-U,,+min{5i,^,}]>-7x/(3r),ie Ji| 

I' AT; (ex) _ >j 

U<^ J2 [Mj{ex) -K{Ai} + mm{Si,6^}]> -jx/{3r),i £ jA 
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u\j2 [BkAliex) + (1 - Bk){Ak{Ex) + U^)] < 7^/(3r)| . 
UeJi ) 

We have the lower bound 

Fj,{A{0,Tf)-CTf>x} 

= Pjo {F^ > x,i e Jo;Gi > x,i e Ji} 

> ¥{Fi >x,ie Jo; Gi >x,ie Ji; -^2(7, 6, 6, e, x)} 

>p|(r-c)A[(£x)- Y. rkNkiex)[n^k} + Sk]> (1 +l)x,i e Jo; 

(r - c)A,(ex) + (r - c)A^,(ex)[E{AJ - ^,] - diiV,(e2;)[E{Ui} + 5i] 

- Y. rkNk{ex)[E{\Jk} + Sk]>{l+7)x,ieJi;E2{7,d,d,e,x) 
keJi\{i} ' 

This probabihty is lower bounded by, for any L [condition on A^j(ex)], 

Y F{E2{i~b~b,e,x)\Ni[ex)=n,,i(iJ^} 

X n ¥{N,{ex)=ni} 

(B.l) xp|(r-c)A[(ex) 

- E rkNk{ex)[E{Vk} + 6k]> {I + 7)x,i £ Ji; 
keJi 

(r - c)Ai(ex) + (r - c)7Vi(ex)[E{Ai} - 5^] 

-d,iV,(ex)[E{Ui} + 5,] 

- E ^fcA^fcM[E{Ufc} + 4]>(l + 7K 
fceJiUi} 

i E Ji I Ni{ex) =ni,ie Ji>. 

Before proceeding, we first state a useful lemma (a proof is given at the end 
of this section). D 

Lemma B.l. For all e, j,5,5 >0, 
(B.2) P{£2(7,^,^,e,^) I Ni{ex) = n^,ie Ji} -^ 1, 
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as X ^> oo uniformly in ni>0,i £ Ji, and 

(^•^) P{ A, > ex} ' ^ 

for all i £ Ji as X — > cx) uniformly in <ni < Lx . 

Equations (B.2) and (B.3) imply that, for any L < oo and ry > 0, one can 
lower bound (B.l) for x large enough by 



ni<Lx,i£j'i K 



- ^ rfcnfe[E{Ufc}+(^fe]n>(l + 7)x,iG Jo; 

(r - c)Ai(ex) + (r - c)ni[E{Ai} - 6i] 
-dini[E{Vi} + 5i] 

- Yl rknk[E{\Jk} + Sk]>{l + -/)x, 

k&Ji\{i} 

i^Ji \Ni{ex)=ni,ieJi } 



X PJ P{Ai > ex}. 

As before, deconditioning upon Aj and A[ and applying a similar change 
of variables as in Appendix A, we obtain the lower bound 

X / l[F\ir-c)A',>il + ^)x+ Yirk-Pk + 6k)yk\ 

X n ^ w - c)^» > (1 + ^)^ + (^* -pi + ~^i)yi 

+ Yl i'^k- Pk + h)yk}dy. 

k€ji\{i} ) 



Now write 

'l<yi<Lx,iGji Jyi>0,i£Ji J{l<yi<Lx,iGji_y 



(1-7?)/ ... = (^i-n)f ...-{i-n)[ 

Jl<Vi<Lx.iGj, Jvi>0.i£j, Jil 
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(the complement taken with respect to the nonnegative orthant). The first 
term in the right-hand side can be handled as in the proof of the upper 
bound (the only difference is the factor 1 + 7 instead of 1 — 7 — (^) . The next 
lemma shows that the second term can be neglected. 

Lemma B.2. 

lim limsup— — ^ / TT P< (r - c)A[ > (1 + 7)3; 

L^oo x^^ P{X) J{l<y,<Lx,ieJiY.^^^^ [ 

+ "^ (i^k - Pk + h)yk 

k£Ji J 



X 



[] pj (r - c)Ai > (1 + 7)x + (d, -pi + k)yi 

+ X! {^k- Pk + ^k)yk \ dy 
keJi\{i} J 



0. 



X 



Proof. The integral over the regions in which at least one yi is smaller 
than 1 is easily shown to be of o(P(x)), so we concentrate on the set {0 < y^ < 

i ^ JiY- The integral 

/ [|pj(^_c)A[>(l + 7)x+ E(rfc-pfc + 4)y4 

n pj (r - c)Ai > (1 + 7)x + {di -pi + 6i)yi 

+ Yl {rk- Pk + 5k)yk \ dy 
kaji\{i] J 

is bounded from above by 
(l[F{{r-c)A:>il+^)x}\ 

V JG Jo / 

E/ l{F\ir-c)A, 



X 

J 



> (1 + 7)x + {di - pi + 5i)yi 
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+ Yl {fk- Pk + Sk)yk, i^Jl> dy. 
k&JAH\ ) 



k&Ji\\i} 

Observing that the integrals can be separated, we obtain the upper bound 

o ( n nK > ^}) E o{n^ > Lx}) n o [ n nK > ^}) 

„ PIA'; > Lx} 

The result then follows immediately. D 

Proof of Lemma B.l. Equation (B.2) follows immediately from the 
following result. Let S„ = Xi + • • • + X„ be a random walk with i.i.d. step 
sizes with E{Xi} < 0. Then 

limsupsupPjSn > x} < lim P< supS„ > x > = 0, 

x^oo „,>1 x^co yn>l J 

since sup„>i S^ is a proper random variable. Apply this result with Xj = 
Ujj — E{Uj} — min{5j, 5j} and Xj = E{Aj} — Aij{ex) — min{^j, (5j}. 
In order to prove (B.3), note that, for rii < Lx, 

FiNi(ex)=ni} ,. ,„ ,. ,r. f o(l)\^'' 

P{Ai > ex} \ X J 

as X ^ oo. The last equality holds because Aj has finite mean. D 
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